Anomalous diffusion through couphng to a fractal environment: 
Microscopic derivation of the "whip— back" effect. 
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Two models for quantum Brownian motion — the Oscillator Bath (OB) model and the Random 
Band-Matrix (RBM) model — are compared and a relation between the spectral density function 
and the variance Vab^ is established. The extension to a fractal environment is then considered 
and the microscopic origin of anomalous diffusion is discussed. In particular, it is shown that 
the asymptotic behavior of the normalized velocity auto-correlation function (VACF) is entirely 
determined by the band form factor. This allows for a microscopic derivation of the "whip-back" 
effect. 
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Anomalous diffusion has attracted considerable atten- 
tion in the last decade (for a review see and nu- 
merous studies have been devoted to its description, 
mostly at the phenomenological level. Among different 
approaches are fractional diffusion equations iQ-^ , Levy 
walks 0, fractional Fokker-Planck equations 



nonextensive statistical mechanics 12 1, quantum Levy 
processes [|l^ and generalized Langevin equations 
Anomalous diffusion often results from the interaction 
with a complex, non-homogeneous background. Exam- 
ples include porous or disordered media [H , chaotic heat 
baths and also fractal environments ||l^ . In this Let- 
ter we focus on a system coupled to a fractal heat bath 
and explore the microscopic origin of anomalous diffusion 
with the help of the RBM model for quantum Brownian 
motion introduced in Ref. ]l7| . In order to make a di- 
rect comparison between the OB model jl^,^ and the 
RBM model, and to extent the latter to the case of a 
fractal environment, we first calculate the bath correla- 
tion function and investigate the validity of the Kubo- 
Martin-Schwinger (KMS) condition. We then relate the 
VACF to the microscopic band form factor by employing 
the first dissipation-fiuctuation theorem and derive the 
"whip-back" effect 

We consider a quantum system S weakly coupled to a 
heat bath B. The coupling is taken linear in the position 
coordinate x of the system. The composite Hamiltonian 
is given by 



H = Hs 



•V 



(1) 



where Hs (Hb) is the system (bath) Hamiltonian and 
V is an operator acting on the bath. It is assumed that 
at t = 0, iS* and B are uncorrelated, so that the initial 
total density operator factorizes, p{0) — ps{0) ^ /5_b(0), 
where ps.sit) — trB,5[/5(t)]. It is further assumed that 
the bath is initially in thermal equilibrium at tempera- 
ture T, pb{0) = ejip{-l3HB), with /3 = (kT)-^. 

In the standard approach to quantum dissipation 
1^,^, the bath is modelled by an ensemble of indepen- 
dent harmonic oscillators (mass m^, frequency LUi) and V 



is assumed to be linear in the positions Xi of the oscilla- 
tors, V = CiXi, where the c^'s are coupling constants. 
The coupling to the heat bath is fully characterized by 
the spectral density function I{uj), which is defined by 



2m,- 



(2) 



On the other hand, in the random-matrix model the 
structure of the bath is not specified, in particular there 
are no environmental oscillators, and ^ is a centered 
Gaussian random band-matrix. In this approach the 
coupling to the bath is uniquely characterized by the sec- 
ond moment Va}? of the random interaction. The explicit 



form used in Ref. |0 reads 



Vah = ^0 {p{ea)p{eh)\ " exp 



2A2 



(3) 



Here e^'s denote the eigenenergies of the bath Hamilto- 
nian (iy^ja) = £a|a)), is the strength of the coupling, 
A the bandwidth and p[e) is the density of states of the 
bath. Using the equivalence of the canonical and micro- 
canonical ensembles in the thermodynamical limit, the 
bath may be parametrized as p7[ 



/Os(0) = \a*){a*\ and p{e*) = p^e 



lie" 



(4) 



where the state |a)* is defined by the temperature T and 
po — Z I \J ixkT'^Cv {Cv being the heat capacity at con- 
stant volume). 

In Ref. we have demonstrated that the averaged 
Markovian master equation for the system obtained from 
the RBM model is identical to the equation derived start- 
ing from the OB model. This shows that the two models, 
although being completely different in nature, are some- 
how related (at least at the level of the master equation) . 
It is our aim to make this connection more precise and to 
directly relate the quantities that define the two models, 
namely the spectral density function (^ and the vari- 
ance (|^). This will be achieved by evaluating the bath 
correlation function. 
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The bath correlation function is defined as 



K{t) = {Vit)V{0))B , 



(5) 



where V{t) — exp{iHBt)V cyip{—iHBt) and = 
trB[pB(0) . . .] denotes the thermal average. Performing 
the average over the random-matrix ensemble in Eq. (|^) 
leads to 



K{t) = 



dsbp{eb)Vab'e'^'- 



(6) 



The above equation shows that K{t) is essentially the 
Fourier transform of the variance Vat^ with respect to 
Eb- On the other hand, the correlation function can also 
be written in terms of the spectral density I{uj) as 

K{t) = J ^^I{uj) ^coth ^^^^ cos(cjt) — ism{ujt) 

^K'{t) + iK"{t) (7) 

by using the identity of K{t) with the kernel of the in- 
fluence functional [|l9|. We see from Eq. that the 
imaginary part K"{t) of the correlation function is sim- 
ply the Fourier sine transform of the spectral density. 
Before we make use of th is ob servation to obtain the re- 
lation between /(w) and Vab^, let us find the most gen- 
eral form of the variance allowed by the laws of thermo- 
dynamics. More specifically, let us determine the con- 
ditions under which K{t) satisfies the KMS condition, 
K{—t) ~ K{t — i/3), which is known to define the ther- 
mal equilibrium state of the bath p^ , ^ . To this end, 
we consider the following form of the variance {k > 0), 



Vab - ^0 [piea)p{eb)] 
Substituting (^) into Eq. i 

Kit) 



exp 



Ac 



P(£a) 



27rA 

2fc^exp 



2A2 
we obtain 



(8) 



(9) 



where we have replaced the density of states p{e) by the 
expression (^. We see from (||) that the KMS condition 
is only fulfilled for the value fc = 1/2 (which is the form 
we have used so far) and that this is independent of the 
actual form of the band. This can be understood by the 
following qualitative argument. Let us write |£a — £b| = 
nd where d = p~^ is the mean level spacing and let us 
assume that nd <^ A. According to (^), the intensity 
of the coupling is approximately given by Aod^''^. As a 
result, wc observe that the coupling is much stronger for 
fc = than for fc = 1/2 (recall that the spectrum of 
B is quasi-continuous which means that d ^ 0). This 
strong coupling then prevents the bath from being in 
equilibrium. In the opposite limit fc = 1, the coupling 
is so weak that the imaginary part K"{t) vanishes. This 
implies that there is no dissipation (see Eq. (p^). 



The most general expression for the variance which 
respects the KMS condition is hence of the form 



Vab = [p{ea)p{eb)] " f(t^Qb = £6 - £a) , 



(10) 



where we have introduced the (symmetric) band form 
factor f(w). Note that the form ( p^ ) ensures the pas- 
sivity of the heat bath in accordance with the Second 
Law of thermodynamics J23] (for fc > 1 the damping co- 
efficient defined in Eq. (nM is negative and the bath is 
clearly active) and also guarantees that the bath remains 
in equilibrium at all times. Now, by comparing the imag- 
inary parts of Eqs. (||) and (0), we arrive at the relation 



I{uj) = 27rsinh ( ^ ) 



(11) 



This equation provides a link between the two models. 
A straightforward application of the above formula using 
the form (||) yields 



I{uo) — 27rAo sinh 











1 exp 


2A2 



(12) 



One may also check that the power spectrum of K'{t) 
satisfies the fluctuation-dissipation theorem |[l^,pj[ 



K'{uj) = coth 



~2~ 



Ku;) 



(13) 



At this point, it is useful to introduce the damping kernel 

il ■ ^^^^ 

In the limit 1 < A < fcT, (|l|) is found to be equal to 

(15) 



M 

With the help of Eq. (|l5|), Eq. ^ thus reduces to 

/(w) - Af7w (16) 

in the limit cj ^ 0. Clearly, the variance (^) defines 
Ohmic damping, as expected from the results of Ref. [0 . 
Let us now turn to the determination of the band form 
factor in the non-Ohmic case. 

In the non-Ohmic regime, the spectral density of the 
bath obeys I{u)) ~ a;",a > 0, at low frequencies. For a 
non-integer, the bath is called fractal |2^]. We immedi- 
ately see that if the band form factor satisfies 

f(c^) ^ |a;r-i , (17) 

in the limit ^ 0, then, according to Eq. (pl|), 

I{uj) - 7r/3 wf(w) - 7r/3 w > . (18) 



2 




1 2 



T3 
C 



1 


(c) 


0.8 




0.6 




/ "A 

/ 0.2 





1 


(d) 


0.8 




0.6 




/ \- 





1 



-2-1 12 
Frequency CO 

FIG. 1. Band form factor f(aj) = I'i'l""^ exp[— c^;^] for four 
values of a: (a) a = 1 (diffusive), (b) a = 2/3 (subdiffusive) , 
(c) a = 3/2 (superdiffusive) , (d) a — 4 (ballistic). 

We conclude that the damping regime is determined by 
the shape of the band at the origin and that the non- 
Ohmic (fractal) regime is the result of a shape of type 
(|l7|). (Notice that the algebraic singularity appearing in 
Eq. (^ for a < 1 can easily be regularized, since the di- 
agonal matrix elements Vaa have no physical meaning, 
by taking Hadamard's finite part x°'~^g{x) dx = 
J^^x"-~^[g{x) ~ g{0)]dx, see e.g. |25|). Furthermore, 
it is known that the coupling to a fractal environment 
gives rise to anomalous diffusion at finite temperature 
In the long-time limit, the mean-square dis- 
placement s{t) — ({x{t) — x{0))^)t for a free damped par- 
ticle grows like for a < 2 and like for a > 2. This 
leads to normal diffusion in the Ohmic case a — 1 and 
to sub(super)diffusion in the case a < 1 (1 < a < 2). 
For a > 2, the growth of s{t) is ballistic. These differ- 
ent regimes can be qualitatively understood by looking 
at the behavior of the band form factor at the origin (see 
Figure]^). Compared to the diffusive case a = 1, we ob- 
serve that subdiffusion is the consequence of the strong 
enhancement of small energy transfers which occurs for 
a < 1 (the absorption by the bath of a given amount of 
energy necessitates, on average, more steps in the latter 
case: the bath becomes "stiff" and diffusion is slowed 
down), while superdiffusion results from the suppression 
of these small transitions for 1 < a < 2. The value 
a = 2 is singled out in all anomalous diffusion phenom- 
ena. This happens here too. When a > 2, the derivative 
at the origin of the function uj°'~^ vanishes and the form 
factor i{Lu) shows a gap around uj — 0. As a result, the 
motion of the free particle is dissipative only as long as 
its energy is larger than the value of the gap. Once its 
energy is smaller, the particle behaves, on average, as if it 
were free (ballistic motion). The foregoing discussion can 
be made more quantitative and the microscopic origin of 



anomalous diffusion more transparent by evaluating the 
VACF of the particle. 

We shall now show that the long time behavior of the 
normalized VACF Cy{t) of the system can be written 
in a simple way in terms of the band form factor. Let 
us recall that the canonical correlation function of the 
velocity operator v is defined as |^ 

{v{0);v{t)) = f]-'^ dX{v{-inX)d{t)) . (19) 
Jo 

In the limit h 0, ( pj| ) reduces to the classical correla- 
tion function {v{Q)v{t)). The first fluctuation-dissipation 
theorem 



CJz] 



(20) 



then relates the Laplace transform Cv[z] of the nor- 
malized VACF Cy{t) = ({i(0);w(t))/(5(0);w(0)) to the 
Laplace transform of the damping kernel j{t), which is 
given by H 



7 z = — 



dui I{uj) z 

We further note that in Eq. (|2l|) 

z 



lim — . 

z^O TT Z'' + LO^ 



= 6{z) 



This gives 



lim l[z] 

z—>0 



t f(-) 
M z^o ^ ' 



(21) 



(22) 



(23) 



where we have used Eq. (p^). Hence, by virtue of the 
final value theorem (assuming that the limit of Ci,(t) 
exists), we can write 



lim a(<) = lim (l + ^z-if(z; 



(24) 



Eqs. (|2^) and ( ]2^ ) show that the band form factor gov- 
erns both the asymptotic behavior of the damping ker- 
nel and of the VACF. To be more precise, we shall as- 
sume that the VACF is of the form Cv{t) ^ for large 
times and employ the following theorem on the asymp- 
totic properties of the Laplace transform of a function 
F{t) II 



If 

then 



ct^, 



Fit) 

zF[z] - cr(A-f l)z-^. 



oo 
■ . 



This theorem is a consequence of 



f^e-^*dt = r(A + l)z-(^+i) 



(25) 



(26) 



In view of Eq. (18), we rewrite the right-hand side of 
©as 
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zC„[z] ~ (l + 







(27) 



Two cases have to be distinguished: a smaher or larger 
than two. In the case a < 2, Eq. (|2^) reduces to 



zCJz] 



MkT 







By using the theorem 



we then infer 



t 



(28) 



(29) 



TT^Q r(a - 1) ' 

For a < 1 (subdifFusion) , the normahzed VACF has a 
negative power-law tail (this can be seen by noting that 
T{a — 1) = r{a)/{a — 1)). This negative correlation 
leads to a incessant change of direction of the velocity 
of the particle and has been named "whip-back" effect 
in Ref. j2^ (see also Ref. for this reason. This ef- 
fect is responsible for the slower diffusion of the parti- 
cle. In contrast, for a > 1 (superdiffusion) , the VACF 
possesses a positive tail which means that the particle 
is more likely to move always in the same direction: 
this results in a enhanced diffusion. In the Ohmic case 
a = 1, the VACF vanishes for large times since r(0) = oo 
and the velocities are therefore uncorrelated. This is of 
course in agreement with the exact form of the VACF, 
Cy{t) = exp(— 7i), which is obtained by taking the in- 
verse transform of Eq. ( po|) with 7[z] — 7. On the other 
hand, in the case a > 2, we note that tends to zero 
in the limit z ^ 0. Eq. (|2^) can therefore be rewritten 
in the form 

zC„[z]-l, z^O, (30) 
and we have accordingly 

Cy{t)r~.l^ t^OO. (31) 

The VACF ( ^l|) corresponds to ballistic motion. 

In conclusion, we have investigated the validity of the 
KMS condition for the RBM model by calculating the 
bath correlation function. We have found that the most 
general form of the variance compatible with the Second 
Law of thermodynamics is given by Eq. (p^. This form 
guarantees the passivity of the heat bath. We have fur- 
ther derived a relation (Eq. (|Tl|)) between the spectral 
density of the OB model and the band form factor and 
extended the latter to a fractal environment. We have 
then discussed the microscopic origin of anomalous diffu- 
sion; first qualitatively by examining the form factor at 
the origin, and second, quantitatively, by expressing the 
asymptotic behavior of the normalized VACF in terms of 
the band form factor. This resulted in the derivation of 
the "whip-back' effect. 

We wish to thank H. Grabert for having suggested the 
direct comparison between the RBM and OB models and 
H.A. Weidenmiiller for useful discussions. 
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